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ABSTRACT 

In  this  paper,  Khater  method  is  used  to  construct  new  analytical  solutions  of  mKdV  equation.  As  a  result,  some 
new  types  of  exact  traveling  wave  solutions  are  obtained  using  trigonometric,  hyperbolic,  exponential  functions  and 
rational  forms.  The  related  results  are  extended  and  the  obtained  results  clearly  indicate  the  reliability  and  efficiency  of 
the  Khater  method. 
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1.  INTRODUCTION 

To  understand  the  nonlinear  physical  phenomena,  arises  in  nature  and  in  the  field  of  applied  sciences  such 
as  fluid  dynamics,  plasma  physics,  solid  state  physics,  optical  fibers,  acoustics,  mechanics,  biology  and 
mathematical  finance,  finding  exact  solution  to  nonlinear  equation  is  important  [1],  In  nonlinear  science,  exact 
solutions  to  nonlinear  partial  differential  equations  play  a  great  role,  since  they  can  provide  much  physical 
information  and  more  insight  of  the  physical  aspects  of  the  problem  and  thus  lead  to  further  applications. 
Wave  phenomena  in  dispersion,  dissipation,  diffusion,  reaction  and  convection  are  very  much  important. 
And  many  kinds  of  effective  methods  have  been  used  to  obtain  explicit  traveling  and  solitary  wave  solutions  of 
nonlinear  evolution  equations  have  been  established,  such  as,  the  ansatz  method  [2],  the  Adomian  decomposition 
method  [3],  the  Darboux  transformation  method  [4],  the  Backlund  transformation  method  [5],  the  inverse 
scattering  transformation  method  [6],  the  Jacobi  elliptic  function  method  [7,8],  the  Exp-function  method  [9,10],  the 
extended  tanh  method  [11],  the  Cole-Hopf  transformation  [12],  the  (G'/G)  -expansion  method  [13-17],  and  the 
modified  simple  equation  method  [18,19]. 

Very  recent  Khater  et  al  [20]  established  a  highly  effective  method  called  the  Khater  method  to  get  exact 
and  solitary  wave  solution  of  nonlinear  partial  differential  equations  (NPDEs.).  This  method  is  among  one  of  those 
general  methods  that  rely  on  auxiliary  equation  and  in  addition,  comprise  seven  methods  namely  (G'/G)- 

expansion  method[21],  improved  (G'/G) -expansion  method  [22],  exp  (-</>(<!;))  method  [23],  extended  tanh- 

function  method  [24],  kudryashov  and  modified  kudryashov  methods  [25],  improved  tan (/^j  -expansion  method 

[26],  and  novel  (G'/G) -expansion  method  [27]  subject  to  certain  conditions  which  makes  it  an  effective  and 
powerful  method.  The  aim  of  this  work  is  to  search  for  new  study  relating  to  the  Khater  method  for  solving  the 
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mKdV  equation  to  demonstrate  the  appropriateness  and  straightforwardness  of  the  method. 

The  organization  of  the  paper  is  as  follows:  in  section  2,  a  brief  description  of  Khater  method  is  given.  In  section 
3,  we  will  study  the  mKdV  equation  by  the  new  methods  and  in  section  4  discussion  on  the  result  are  provided.  Finally 
conclusions  are  given  in  section  5followed  by  references  in  section  “References”. 

2.  SUMMARY  OF  KHATER  METHOD 

In  2017,  Khater  method  was  first  proposed  by  Khater  et  al  [20]  for  solving  nonlinear  partial  differential  equations. 
Firstly,  we  consider  a  general  nonlinear  evaluation  equation  in  two  independent  variables,  say  X  and  t , 

(p(u,ux,utuxxuttuxt, . )  =  0,  (2.1) 

Where  U  —  u(x,  t)  is  an  unknown  function  to  be  determined,  tp  is  a  polynomial  of  Li(x,t  )  and  its  partial 
derivatives  wherein  the  highest  order  partial  derivatives  and  the  nonlinear  terms  are  involved  and  the  subscripts  indicate  the 
partial  derivatives.  The  key  steps  of  the  method  are  as  follows: 

Step  1:  We  now  introduce  the  traveling  wave  variable, 

£  =  X-Vt ,  U  -  u(x,t)  =  «(£)  ,  (2.2) 

where  V  is  the  speed  of  traveling  wave  and  the  wave  variable  (2.2)  transforms  Equation  (2.1)  into  the  following 
ordinary  differential  equation  (ODE): 

(j)(u,U  ,U  ,U* , . )  =  0,  (2.3) 

where  (j)  is  a  polynomial  of  ll  and  its  derivatives  and  the  super  scripts  stipulate  the  ordinary  derivatives  with 
respect  to  £ . 

Step  2:  In  many  instances.  Equation  (2.3)  can  be  integrated  term  by  term  one  or  more  times,  yielding  constants  of 
integration,  which  can  be  set  equal  to  zero  for  straight  forwardness. 

Step  3:  The  Exact  solution  of  Eq.  (2.3)  is  supposed  as 

i v 

u(%)  =  Yjaiaf(S) ,  (2.4) 

i= 0 

where  ar  a  ,  is  constant  to  be  determined,  such  that  tlN  ^  0  and  f  (£)  satisfies  the  following  second  order  linear 
ordinary  differential  equation  (LODE): 

/'(£>  =  (aaf(i>  +jB+ <ram ) ,  (2.5) 

ln(a) 

Step  4:  The  positive  integer  N  can  be  determined  by  using  the  homogeneous  balance  between  the  highest  order 
derivatives  and  the  nonlinear  terms  appearing  in  (2.3). 

Step  5:  Using  Eq.  (2.4)  and  Eq.  (2.5)  into  Eq.  (2.3),  collecting  all  the  terms  having  same  powers  of  , 
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where  (z  =  0,  1,  2,  . )  and  equating  to  zero,  weform  a  system  of  algebraic  equations.  This  system  is  further 

solved  symbolically,  to  determine  the  values  of  a'  ,CX ,  /?  and  (T . 

The  solutions  of  Eq.  (2.5)  subject  to  the  coupled  cases  are  given  as  follows 


If  (J3  -aa<0&a*0). 


Jit) 


P  J(P  -  aa) 

-  —  + - tan 

a  a 


-iP  -aa) 


t 


or 


jit)  _ 


P  J{P  -  aa) 

-  —  + - cot 


-(p  -aa) 


If  (Pz  -  aa  >  0  &  a  ±  0). 


Jit)  _ 


P  JP  -aa) 


tanh 


-(p  -aa) 


t 


Jit) 


P  JP  -aa) 


coth 


-(P  -aa) 


t 


If  (P~  +aZ>0&a^0&a  =  -a). 


Jit)  _ 


P  Jp~  +  a2) 

- h  — - tanh 

a  a 


iP2  +  a2) 


Jit)  _ 


P  Jj2  +  a2) 

— I - coth 


a 


a 


02  .  2  ^ 


iP  +a  ) 


t 


If  {P2  +a2  <0&a*0&a  =  -a). 


Jit)  _ 


P  J{p2+a2)t 

+  - - tan 


a 


a 


V 


-(P2  +  a2) 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


(2.11) 


(2.12) 
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(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 

(2.20) 
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flit  _ 


If  (or  =  a  =  0). 


Jit)  _ 


-(1  +  e 


2^)±^ 


+  1) 


jit) 


-(1 +e^y. 


2^)±i4^+6^+1 


2e 


m 


If  (J3  =  aa). 


Jit) 


-a(P$  +  2) 

P2Z 


If  (fi  =  k,a  =  2k,  a  =  0). 


Jit) 


e  r  -1 


If  ( fi  =  k ,  cr  =  2k ,  a  =  0). 

jit)  r 


l-e 


w 


If  (2/3  =  a  +  a). 


Jit) 


1-ore2 


Ua-ctJ 


I 

l-<7e2 


(a-aj 


Jit)  _ 


^(ar-a)# 
orez  +1 


1 


-ere2 


(or-er)£ 


-1 


If  (-2/3  =  a  +  a). 


(2.21) 


(2.22) 


(2.23) 


(2.24) 


(2.25) 


(2.26) 


(2.27) 


(2.28) 
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a 


m 


\(a-o)£, 


+  a 


+  a 


(2.29) 


lf(«r  =  0). 


J  (£) 


fie 


M 


1+ °-eK 
2 


If  (p  =  a  =  a  ^  0). 


am 


-(a%  +  2) 
at, 


If  (p  =  cr  =  0). 


af®  = 


If  (p  =  a  =  0). 


a 


m 


If  (fi  =  0,a  =  <j). 


am 


tan 


a^  +  C  j 


If  ((7  =  0). 


am 


a 

ip 


(2.30) 


(2.31) 


(2.32) 


(2.33) 


(2.34) 


(2.35) 


Where  C  is  an  arbitrary  constant. 

Step  6:  Substituting  the  above  values  and  the  solutions  of  Eq.(2.5)  into  Eq.  (2.4)  weattain  the  exact  solutions  of 
Eq.  (2.1). 

3.  KHATER  METHOD  IN  mKdV  EQUATION 

The  mKdV  equation  arises  in  nonlinear  optics  and  in  the  propagation  of  long  internal  waves  in  a  fluid  when  the 

2 

coefficient  of  the  usual  nonlinear  term  in  the  KdV  equation,  uux  is  zero  and  the  higher  order  nonlinear  termw  ux  dominates 
over  higher  order  dispersive  terms.  We  consider  the  mKdV  equation  in  the  form: 

UT  +  Cu2Ux  +  puxxx  =  0 ,  (3.1) 
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Where  u  —  and  £  ,  p  are  non  zero  constants.. 

Now,  using  the  wave  variable  u(x,t )  =  u(£) ,  £  =  X  —  Vt ,  in  (3.1)  and  integrating  the  resulting  equation  and 
neglecting  the  constant  of  integration  we  find 


-Vu  +  -£u3  +  pu"  =  0, 

3 


(3.2) 


By  balancing  U3  and  U  we  get  3./V  =  iV  +  2 ,  henceiV=l.  We  then  suppose  the  Eq.  (3.2)  has  the  following 
formal  solution 


u{£)  =  a0  +  alaf^) , 


(3.3) 


/(<f) 


Substituting  Eq.  (3.3)  into  Eq.  (3.2),  we  use  (2.5).  The  left-hand  side  of  (3.2)  becomes  a  polynomial  in  a 
Setting  the  each  coefficient  of  different  power  of  a‘^h)  ,{i  —  0,1, 2, 3 . )  to  zeroyields  a  system  of  algebraic  equations. 


£ Aq  -  3vaQ  +  3 ajfpa^  =  0 


3  2 

If  a  q  +  6pa  Oj  =  0 


-  Vflj  +  lapociy  =  0 

2 

faQ  a£  +3j8p<7ciy  =  0 

Solving  the  algebraice  quations  by  Maple  we  get 


1  yj-6£pfi  sl~6Cpo-  1  2, 

- ,  a,  = - ,  v  =  -  p(4aa  -  B  ) 

2  tf  X  £  2 


So  that,  the  exact  traveling  wave  solution  of  Eq.(3.2)  is  of  the  form: 


2  c  t 

which  leads  to  the  following  travelling  wave  solutions  subjected  to  the  corresponding  cases 


(3.4) 


(3.5) 


If  (jB  -act  <0&(T  *0) 


u(x,  t)  = - 

2  £ 


1  sl-6fpj3  J-6 fp  [  _  n2 
T - r - +  — - — V'  ' 


aa-  tan 


i 


aa-  [3 


■(x-vt) 


y 


(3.6) 


or 
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iftCpp  fefr  I  ~2 

u(x,t)  = - 1- - vena-  B 

2  C  C 


cot 


a<r-  p~ 


(x-vt) 


J 


If  (J32  -aa>0&a7= 0). 

,  n  1  f^pP 

u(x,t)  =  -- 


2  c  C 


-  cm  tanh 


i? 


-aa 


{x-vt) 


J 


u(x,t )  =  - 


1  f6frp 


2  <r  <r 


—  acr  coth 


-aa 


(x-vt) 


J 


If  (p2  +  a2  >0&ar=0&a  =  -a). 


1  4-6&P  yl-SCpc 

u(x,t)  =  — 


2  c  C 


P  i?2  +«2  , 

— I - tanh 

a  a 


\ 


(x-vt) 


J 


or 


u(x,t)  = 


1  yl-6frfi  sl~6Cpa 


2  C  C 


P  J?2+a2 

—  + - coth 

a  a 


(x-vt) 


J 


If  (P2  +  a2  <0&<7^0&£T  =  -a). 


1  V-6 TpP  ?6Cpa 

u(x,  t)  =  — 


2  {  { 


a  a 


2  2  ) 
+  a 


-  tan 


2  2 1 
+  a 


-(x-vt) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


or 


u(x,t) 


1  yJ-6CpP  |  V~6  jpo_ 

2  c  C 


(3.13) 


If  (P2 -a2  <  0 &  a  =  a). 


u(x,t) 


1  V~6  CPP  |  V~6  CP° 

2  C  C 


V 


J 


(3.14) 
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or 


1  V-6  CpP  V-6 ~Cpct 

u(x,t)  =  — 


2  c  C 


- - cot  - - 


a  a 


-( x-vt ) 


(3.15) 


If  (/?~  -  aT  >  0  &  a  =  a). 


,  A  1  -l-BCpfi  -j-BCpcr 
u(x,t)  = - 1- 


2  C  i 


J-+ 


(f-a2) 


tanh 


a  a 


(f-a2) 


(x  -  vt) 


(3.16) 


or 


,  „  1  s]-6Cpfi  ^ 

u(x,t)  = - 1 - 

2  (  i 


coth 


a  a 


(f-a2) 


(x-vt) 


(3.17) 


If  (act  <0&£r^0&/?  =  0). 


u(x ,  t )  = 


\l-6Cpc 


tanh 


au 


\ 


-(x-vt) 


7 


(3.18) 


u(x,t )  = 


V-6 Jpo 


coth 


T- 


aa 


\ 


(x-vt) 


7 


(3.19) 


If  (fi  =  0  &  a  =  -a). 


u(x,  t )  = 


u(x,  t)  =  ■ 


V-6  CP<t 

-(!  +  e2flr(x-V,) )  ±  ^2(e4a(x-Vt) +l) 

C 

2  a(x-vt)  , 
e  - 1 

sl~6Cp<y 

2 a(x-vt)..  1  4 a(x-vt)  ,  2a(x-vt)  . 
-(1  +  e  )±\e  +6e  v  '+1) 

C 

„  2 a(x-vt) 

2e 

(3.20) 


(3.21) 


If  (P2  =  aa). 


s]-6Cp/3  ^i-eCpcr 

-a(p(x-vt)  +  2) 

c  c 

p2(x  —  vt) 

(3.22) 
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If  (fi  =  k,a  =  2k,a  =  0). 


Ji(x-vt) 


1  k(x—vt) 


If  (2 J3  =  a  +  a). 


.  1  V-6 ~CpP  sl~6Cpv 

u(x,t)  = - 1 - 

2  c  c 


\(a-o)(x-vt) 
1  -aeA 


\(a-<j)(x-vt) 

\-ael 


1  yj-6CpP 

u(x,  t)  = - 1 - 

2  c  C 


,  1  V-6 ~CpP  sl-6Cpo- 

u(x,t)  = - + 


2  C 


C 


Ua-o)(x-vt) 
ae z  +1 


(a-a)(x-vt) 


-ae" 


-1 


If  (~2fi  =  a  +  a). 


,  A  1  Jpfi  ^-6Cpa 

u(x,t)  = - 1 - 

2  c  c 


1 


(a-a)(x-vt) 


+  a 


Ua-o)(x-vt) 
ez  +a 


If  (a  =  0). 


,  A  1  /-6 TpP  -J-bCpv 

u{x,t)  = - 1 - 

2  c  c 


fieP {x~Vt) 

1+"efi<x~vt) 


If  (p  =  a  =  a  *  0). 


1  V-6  frP  ^/=6^^J 

u(x,  t)  = - - 

2  c 

If  (j3  =  a  =  0). 


C 


-(a(x-vt)  +  2) 
a(x-vt) 


u(x,t)  = 


J-bCpe 


-2 

a(x  -  vt) 


If  (P  =  Q,a  =  a). 


u{x,t)  = 


tan 


r  a(x-vt)  +  C^ 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


(3.30) 
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4.  DISCUSSIONS  ON  RESULT 

Different  types  of  solutions  of  the  mKdV  equation  obtained  by  Chen  et  al  [28]  using  double  sub-equation  method, 
Triki  and  Wazwaz  [29]  applying  sub-ode  method,  Wanget  al  [30]  implementing  the  (C/  /  G)  -expansion  method,  Yan  [31] 

utilizing  the  extended-tanh  function  method,  Smyth  and  Worthy  [32]employing  the  approximate  method,  Biswas  et  al 
[33]using  Lie  symmetry  method.  But  implementing  the  suggested  Khater  methoda  huge  amount  of  travelling  wave 
solutions  of  mKdV  equation  have  been  generated,  which  consist  of  trigonometric  function,  hyperbolic  function, 
exponential  function  and  rational  forms.  We  realize  that  our  established  solutions  of  mKdV  equation  contains  the  solutions 
found  by  the  aforesaid  methods  as  well  as  different  new  solutions,  which  obvious  that  Khater  method  is  more  general  one. 

5.  CONCLUSIONS 

We  have  studied  the  mKdV  equation  using  Khater  method.  And  established  a  rich  class  of  general  new  exact 
travelling  solutions,  comparing  the  existing  solutions.  The  new  found  solution  might  be  significant  to  analyze  the 
complicated  phenomena  arising  in  the  areas  of  physics,  engineering  and  other  physical  sciences  and  able  to  expand  the 
influence  on  future  research  work.  Moreover,  our  study  shows  that  Khater  method  is  effective,  consistent,  straightforward 
and  practically  well  suited  in  finding  exact  traveling  solutions  as  well  as  different  types  of  partial  differential  equations. 
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